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The quark hadron phase transition in a rotating compact star has been studied. The NLZ-model 

has been used for the hadronic sector and the MIT Bag model has been used for the quark sector. 

It has been found that rotating twin star (third family) solutions are obtained upto ft « 4000s~^. 

^h ' Stars which are rotating faster than this limit do not show any twin star solution. A backbending 

(«_^ . in moment of inertia is also observed in the supermassive rest mass sequences. The braking index 

^~; I is found to diverge for a star having pure quark core. 

a- 

I ^ ■ One of the best possible laboratories to study strongly interacting matter at super-nuclear densities is the compact 

_^ ,: stars [1]. The matter density near the core of a compact star can be even 10 times that of normal nuclear matter. 
fSj At such a high density different exotic phase transitions may take place in the strongly interacting matter [2]. Some 

of them are the quark-hadron phase transition, the kaon condensation and hyperonic phase transition [3]. There has 
^ . been a large number of efforts to correlate these phase transitions with the observable properties of a compact star. 
^\ ] However, all the stars, observed so far, can be explained within the periphery of strongly interacting nuclear matter 

without invoking any exotic phase transition. 



o 
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\^^ ■ For a static (non-rotating) compact star the one possible signature is to look at the phase transition is the mass- 

CD , radius relationship. However, the situation is little bit different for a rotating star. It has been argued by Glendenning 
et. al. [4] that one of the important quantity to probe, in order to study the phase transition, is the braking index. The 
^ central density of a rapidly rotating pulsar increases with time as it spins down and when a certain central density 
Q,' is achieved a quark phase may appear at the core of the star. Glendenning et. al. [4] found that when the pure 
J,^ I quark core appears the star undergoes a brief era of spin-up and the braking index shows an anomalous behaviour. 
?H ■ Furthermore, when the moment of inertia was studied as a function of angular velocity it showed a "backbending". 
c/3 ■ Motivated by the above results some authors have studied the backbending phenomenon using different EOS and 

• • ,: also with different rotating star codes. Recently Spyrou and Stergioulas [5] have published an interesting result. They 
. !_^ I have studied the rest mass sequences for a particular EOS (with a quark-hadron phase transition) using their code 
^^ . "rns". They argued that previous results needed some numerical refinements and concluded that backbending is 
J3 i observed only for the supermassive sequences. 

Recently some of us have used the NLZ model (which is a variant of the non-linear Walecka model) for the hadronic 
sector and MIT Bag model for the quark sector to look at the quark-hadron phase transition in static compact stars 
[6] . In that work it was found that there was a solution for the third family of stars known as twin stars [7] . In this 
article we employ the same EOS for a rapidly rotating star. The basic motivation is to study the fate of the twin 
stars in a rotating model and also to study the possible signatures of these stars. First we will describe the model 
that we use here. Then the General Relativistic features of a rotating star will be briefly outlined and at the end we 
will discuss the results. 

In this paper we will use one variant of the non-linear Walecka model, called the NLZ model [8], for the hadronic 
sector and the MIT Bag model for the quark sector. 

The lagrangian density for the hadronic model that we use is given by [8] : 

£^Co+Cyy+Ci (1) 

where 



-Co = ^ V's (*/9- me) ^B + \d*'od^,G ^ U (a) - ^G^'^'G^, + U {to) 

B 

- Ib^-'B^,, + \mlR^R^ -J2^B (9<rB<J + 3c.s^^7m + 3p^^7p^b) V'b (2) 
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Cyy ^ 1 {d'^a*d^a* - ^ml.a*') - ^^''S^, + ^m^^^^^ 

- X! ^S {ga'Bd* + g4,B<l>''7l^) ^B (3) 

B 

A= ^ V^(i;9-mO^/ (4) 

/— e,/i 

In the above equations the ^^ runs over all the baryons (p, n, A, E", S+, S^, S° and S^) and the ^^ runs over 
all the leptons (electrons and muons). The piece of the Lagrangian Cyy is responsible for the hypcron-hypcron 
interactions [8]. The meson fields are cr, w, i?, (T*(/o(975)), and <p. The sigma and omega meson potentials are given 
by [8-10] 

_ {m^of .92 g' 33 fT* ,.. 

U. - ^^— + ^- + ^ ' (5) 

^"" 2 4 ■ ^^ 

The nucleon coupling constants are chosen from the fit to finite nuclei properties. The vector coupling constants of 
the hyperons are chosen according to the SU(6) symmetry as [11-14] 

1 1 1 

-zQi^N = -^gujA = xfft^S — 5wS 
1 

9pN - -zgpT. — gps 
gpA = 

2V2 

2g</.A — 2g^s — g<t>B. — — 77-gi^N 

g<pN = (7) 

The hyperonic scalar coupling constants are chosen to reproduce the measured values of the optical potentials [15,16], 

U[^^ ^ U'^^ = -30Mey, t/|^^ = -28Mey 

\jf^ = \jf^ = 2v'f^ = 2vf^ = -40Mey (8) 

We use the parameter sets as given in ref. [6] 

At the mean-field level, the meson fields are replaced by their ground state expectation values. The equations of 
motion for different meson fields then can be obtained by standard methods, they are given by 

d 



mlcr + -p-U{(T) ^ ^ g^Bpl 
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^l'^* ^^g^'Epf 



mlujo + gzu^i^l = ^ 5c^B 



mlRzfi = ^gpBT: 
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micjya 



Y.d'^BTip^ (9) 



The single particle energies for baryons follow from the Dirac equation 

Esik) = gsiu^o + gB4,4>o + gspT^ po + y /c^ + m*^ (10) 

where the baryon effective masses are given by 

m*g = mB + gaBcr + ga'Bcr* (11) 



For Icptons the energy is Eik) = \Jk'^ + mf^ and we work with the vacuum masses of the leptons. 
The pressure and energy density obtained from these models can be written as [6] 

1 2 2 I .92 3 , .93 4 I 1 2 *2 ,1 2 2,3 . 



\mlRl^ + \ni\<g + Y. (^ / " ^^^ 



k^ + mf (12) 
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k^+mf 
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where Vi is the degeneracy factor of the i-th species. This EOS will be used to look at the star properties. 

For the quark sector we use the MIT Bag model [17]. The pressure and energy density of the quark matter is given 

by 



Y. i^ I dkkymj + k^ + B (14) 

f—u,d,s 



gO - Y^ ^ / rlhh^.f^ 



P'^= y. ^ dkk' , -B (15) 

In the above two equations the sum runs over the three flavours of quark and Vf is the degeneracy. 

In this work the deconfinement phase transition is assumed to be of first order which proceeds via mixed phase. At 
zero temperature the mixed phase should follow Gibbs criterion [4] in presence of two conserved charges. The Gibbs 
criterion dictates 

P"{^lb.|Ie)=P'^{^h.^^e) 

Mb - Mf - A^? (16) 

A^e = Mf = A^? (17) 

where P^, ji^ and ji^ are the hadronic pressure, hadronic contribution to the baryon chemical potential and hadronic 
contribution to the charge chemical potential respectively. The similar quantities for the quark phase are denoted as 
P'5, ^^ and pf^ for the quark phase. 

The volume averaged energy density in the MP can be written as 

e-(l-A)e^(/i6,Aie) + Ae^(Ai,„Aie) (18) 

where A is fraction of quark matter present in the mixed phase. 

In the quark sector, we have taken light quark masses to be zero and strange quark mass to be 150 MeV. So bag 
pressure B is the only parameter. The hadron-quark phase transition is possible for B^^^ in the range 175 - 190 MeV. 
But twin star solution is obtained for a narrow range for B^/^ « 180 - 182 MeV. In the present letter we have given 
the resuhs for B^/^ = 180 MeV. 



The EOS is plotted in figure 1. From the figure one can see that the phase transition starts at around ec — 
4 X Vd^^gmjcrn? and it ends around ec = 1 x VS^^ gvn j cm? . We would like to mention at this stage that due to the 
phase transition there is a substantial change in the slope (sound speed) of the eos. This point will be important in 
the context of our results discussed later. 

Once the EOS is obtained the next job is to solve the Einstein's equations for the rotating stars using the EOS. To 
solve the Einstein's equations we follow the procedure adopted by Komatsu et.al. [18]. In this work we briefly outline 
some of the steps only. The metric for a stationarily rotating star can be written as [19] 

ds^ = -ef+Pdt^ + e^" (rfr^ + r^dO^) + e'^-Pr^sin^9 {d(j) - udtf (19) 

where a, 7, p and lu are the gravitational potentials which depend on r and 9 only. The Einstein's equations for 
the three potentials 7, p and w have been solved by Komatsu et.al. using Green's function technique. The fourth 
potential a has been determined from other potentials. All the physical quantities may then be determined from 
these potentials [19]. 

Solution of the potentials, and hence the calculation of physical quantities, is numerically quite an involved process. 
There are several numerical codes in the community for this purpose. In this work we use the 'rns' code. This code 
developed by Stergioulas is a very efflcient in calculating the rotating star properties. We do not mention the details 
of the code here; they may be obtained in ref. [18,19]. We have used this code to study the rest mass sequences and 
the O sequences for both normal and supermassive stars. 

In figure 2 we show the mass-radius relationship for different values of Q starting from the static limit to f7 = 
5500s~^. From the plots one can see that as the value of fl increases both the radii and the masses increase, which 
is an obvious result. However we find some interesting results. The static mass radius plot shows that there are 
occurrences of twin stars in this model. The same interesting result, i.e. the occurrence of twin stars, is also observed 
for rotating stars. However, as we increase n, the third family solution becomes less probable as can be seen from the 
plots. The third family solutions are obtained until the rotational velocity is « 4000s^^. For stars rotating at higher 
angular velocity, the third family of stars cannot be observed. 

In figure 3 we have plotted the moment of inertia (/) as a function of Q for different rest masses i.e. for both 
normal and supermassive sequences. For the normal sequence the moment of inertia increases with the angular 
velocity monotonically. However it is not exactly the same for the supermassive sequences. For a supermassive 
sequence, there is one branch of the curve for which the moment of inertia increases with angular velocity but there 
is another branch for which the moment of inertia decreases with increase in f2. This phenomenon is known as 
backbending. As pointed out by several authors, this phenomenon could be a possible signature of the quark-hadron 
phase transition. Here we would like to emphasize that even for supermassive sequence, the backbending is only 
observed for the cases where density in the core is high enough to facilitate pure quark phase. For example though 
the range of rest mass 1.35Mq < Mq < 1.48Mq lies in the supermassive domain it does not have a pure quark core 
and it also does not show any back bending. 

In the next figure i. e. in figure 4 the angular momentum is plotted as a function of energy density for different rest 
masses. We have plotted these curves only for the supermassive sequences. From the plots we see that the angular 
momentum (J) decreases with ec then it increases (which is the unstable region) and then again decreases showing a 
third family solution. 

In figure 5 the braking index has been plotted against Q. As pointed out by Glendennig et al the divergence of 
braking index is a signature of the phase transition. This quantity is defined as 

The braking index is found to diverge for supermassive sequences whereas no fluctuation has been observed for 
the normal sequence. The amount of the fluctuations is a manifestation of the size of the quark core of the star. 
For a star in the normal sequence, the central density is such that the core of the star lies at the mixed phase and 
hence the braking index does not show any fluctuation. For the supermassive sequence, in the stars for which the 



central density is high enough to accommodate a quark core, the braking index diverges. For the normal sequence 
that we have plotted here the maximum central density is ^ 8 x lO^^gm/cnV^ whereas for the supermassive sequence 
it is ~ 2 X lO^^gm/cm^. Comparing these numbers with the different central energy densities as shown in figure 1 
it is quite obvious that for the normal sequence the core is in the mixed phase region whereas for the supermassive 
sequence the core of the star consists of pure quark matter. 

To summarise, we have studied the occurence of twin solutions in rotating compact stars. The third family solution 
or twin stars occur due to a substantial change in the sound speed. We have found that though the twin solutions 
are present in rotating stars they vanish after a certain value of fl. Hence stars which are rotating with very high 
angular velocity do not show a twin solution. There is another interesting observation in this work compared to that 
of Spyrou et.al. [5]. They have found that there is a small oscillation of the braking index in the normal sequence. 
However we do not find any such oscillation in the normal sequence. This may be due to the fact that in the central 
region we have only a mixed phase whereas they have a quark core in the central region. This is obviously ascribable 
to the difference in the EOS. 
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FIG. 1. Equation of state of the model considered here. 
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FIG. 2. Radius-Mass plots for diflferent n. 
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FIG. 3. Moment of inertia as a function of n for different rest mass. 
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FIG. 4. J as a function of tc for different rest mass. 
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FIG. 5. Braking index as a function of fl. 
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